Is the tetraneutron a bound dineutron-dineutron molecule? 
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Abstract 

In light of a new experiment which claims a positive identification, we discuss the possible 
existence of the tetraneutron. We explore a model based on a dineutron-dineutron molecule. We 
show that this model is not able to explain the tetraneutron as a bound state, in agreement with 
other theoretical models already discussed in the literature. 
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In a recent experiment Q| , the existence of bound neutron clusters was studied by frag- 
mentation of intermediate energy (30-50 MeV/nucleon) ^^Be nuclei. In particular, the frag- 
mentation channel ^°Be + ^n was observed and the ^n system was tentatively described as 
a bound tetraneutron system. 

The possible existence of the tetraneutron has been discussed theoretically by numerous 



authors already in the 1960's (e.g., refs. 



that time 
with ^He 
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). The experimental search at 



including the experiments on double pion charge exchange reactions 



1J,|15|, gave negative results. The general conclusion, based on experiment and 



theory, was that the tetraneutron cannot be bound. More recent experiments |l 
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21 



3, 



13, 



24| give support to this assertion, except perhaps 



and theoretical calculations 

ref. 2^ where a hyperspherical function method was used and led the authors to the 
conclusion that the tetraneutron may exist as a resonance in the four-body continuum at 
energy of about 1-3 MeV. These theoretical models have recently been reconsidered in ref. 
j3| where it was concluded that a too strong four-nucleon force is needed to bind the 
tetraneutron and that this force would unreasonably bind ^He by about 100 MeV. 

It thus seems very unlikely that the tetraneutron can be explained within any standard 
theoretical model. In this brief report we consider the possibility that the tetraneutron can 
be described as a composite dineutron-dineutron molecular system. To our knowledge this 
is a hitherto unexplored model for the possible binding of neutron matter. In this model, 
the dineutron is considered to be slightly bound due to a polarization mechanism induced by 
the presence of the other dineutron. Although the dineutron system is not bound (it has a 
virtual singlet state at = 66 keV), there is some theoretical (221, |2Sl, |2^ and experimental 
[30I 1 evidence that it might become bound in the presence of another nuclear system. 

Another reason to believe that neutrons clusterize, or form correlated bound pairs inside a 
nucleus, is the existence of the Borromean nuclear systems, e.g. ^He and ^^Li, or nultineutron 
configuration systems in light nuclei, e.g. ^He. The Borromean systems are thee-body 
systems in which the particles are not bound in pairs, but the three-body system is bound. 
For example, although ^''Li and the two-neutron systems are not bound, ^^Li is bound. It is 
the presence of ^°Li which induces the neutron-neutron binding correlation, as first suggested 



by Migdal [27| . The authors of ref. [321] also suggest that two of the three most likely ground 
state configurations in ®He are compatible with an a+'^n cluster system. 

The model for the tetraneutron is based on two clusters of dineutron molecules as dis- 
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played in Fig. A similar model was used in the framework of the generator coordinate 
method for the scattering problem in ref. ISSll. The molecular models, similar to those 
used in quantum chemistry, were successfully applied to the chain of Be isotopes j3j]. The 
dineutron wave function can be written as 



(y9i2 = ^Avi \ No exp 
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(1) 



where = 2 |ri + R/2| is the (relative) intrinsic coordinate and ri is the position of the 
neutron 1 with respect to the center of mass of the tetraneutron. The relative motion 
spatial wave function is taken as a Gaussian with A^o = (fr^Tr^/^^ ^ ^^d b is the oscillator 
parameter. The dineutron is assumed to be in a spin-0 state, described by the spinors Xi^^ 
(spin-up) and xi (spin-down), respectively. The wave function is antisymmetrized by the 
operator Au = | (1 ~ -^12)? where P12 exchanges neutrons 1 and 2. An analogous wave 
function, ^934, is written for the second dineutron molecule. Since the spatial part of the 
dineutron internal wave functions is even under inversion, the operator Au only acts over 
the spin variables and we can write Eq. (|T)) as 
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^{Xl X2 



AO) 
A12 • 



(fi2 = Nq exp 
= A'o exp 

Accordingly, the wave function for the second dineutron is 

9234 = A"o exp 
The total wave function of the tetraneutron is 



(2) 
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(0) 
X34- 



(3) 



^ = VGAuM [<^12¥'34]$(R), 



(4) 



where the operator ^12,34 = | (1 — P13 — Pu — P23 — -P24 + -P13-P24 + -PL4-P23) implements the 
antisymmetrization between the dineutrons. The factors in Eq. and in Eq. 
(jl} account for the proper normalization. The function $ (R) is the dineutron-dineutron 
molecular wave function. 

Using Eqs. (1-4), we obtain (here we drop the upper index (0) in the singlet spinors) 



^ = Vgn^ {/ (^g, ^b) X12XU - g (^g, ^b, H.) X14X32 - h (C„ Cfe, R) X13X24} $(R), (5) 
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FIG. 1: Tetraneutron as a dineutron-dineutron molecule. 



where 



and 



9 {iai ^b, R-) = 3/ i^a, ^b) exp 



h (Ca> = 3/ (^a, ^fc) exp 
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The total Hamiltonian for the tetraneutron system is 



H 



2m 



N 



i=l 



The neutron-neutron potential was taken as a two-body Volkov potential 



(6) 



(7) 



(8) 



(9) 



V = J2{^- M + MP^^ Vi 



V,, (r) = exp (-rVa^) + Vf, exp (-rV/?^) , (10) 



where Pfj is the Majorana exchange operator and the parameters Va, Vp, a and (3 are 
chosen to reproduce the scattering length and effective range in low energy nucleon-nucleon 
collisions, as well as the binding energy of '^He 
The total Hamiltonian H can be presented as 



H = Ta + n + Tn + V, 



(11) 



with 



Ta,b — --. At , Tr — ^j^. 



-A, 



(12) 




FIG. 2: The effective potentials Ui{R) and U2{R) entering Eq. 1)20(1 . The sohd curve is the sum of 
the two potentials. 

Since in our model the mean positions of the clusters are taken to be at R/2 and — R/2, 
respectively, and since the position of each nucleon is symmetric with respect to the origin 
of each cluster, there is no admixture of spurious motion of the center of mass. 

We use Eq. (jSJ as a starting point of a variational procedure. To obtain an effective 
Schrodinger equation for $ (R), we multiply the full 4-body Schrodinger equation if \1/ = E"^ 
by "^y^ (R) and integrate over and We will assume the variational function $ (R) 
to be independent of angles (i.e., an s-wave state), 

*(R) = ^^fl. (13) 

v47r 

so that the operator Tr in eq. (fT^ becomes 

The integrals over and ^2 can be performed analytically. One obtains an effective 
Schrodinger equation for u{R) in the form 

u"{R) + RV3{R)u'{R) + {Vi{R) + V2{R) - E} u{R) = (15) 

2?77,jv 
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FIG. 3: (a) The total potential Ui{R) + U2{R) entering Eq. (|20|) for 8 different parameterizations 
of the Volkov potential. The oscillator parameter b = 1.5 fm was used, (b) The same as in (a), 
but using the Volkov-1 interaction and varying b from 1.2 fm to 2.0 fm. 
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Since we are interested only in the relative energy between the clusters, we extract from 
the total Hamiltonian the internal kinetic energy of each cluster: T*"* = ?)%^ / {Amb'^) and 



V'''' = 14 (1 + &V"^) + Vp{l+ 62//32) that do not depend on R. 
Eq. can be put in a conventional form of a Schrodinger equation: 

^2 



3/2 



2mAr 



t;"(i?) + {f/i(i?) + U2{R) - E}v{R) = 0, 



(20) 



where 

v{R) = u{R) exp 
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f{R')dR' 

The effective potentials in eq. 1201 are given by 

U2{R) = V2{R), (22) 

and 



(23) 



We use for the oscillator parameter b = 1.5 fm. A set of 8 parameters of the Volkov 
potential were taken from Ref. j^. In Fig. |21we show the results for the Volkov VI force. 
We notice that the potential U2{R) is repulsive and the potential Ui{R) is attractive. Their 
sum is dominated by the repulsive part. 



In Fig. El^a) we show the total potential 
the Volkov potential, as taken from Ref. 



Ui{R) + U2{R) for the 8 sets of parameters in 



36|. One sees that none of the parameter sets 



leads to a potential with an attractive pocket, which could be a sign for the existence of a 
bound state. But even in that case the pocket would have to be deep enough to allow for 
the appearance of the bound state. 

Finally, we have varied the oscillator parameter from 1.2 fm to 2 fm for each set of 
parameters of the Volkov interaction. No pocket appeared in the effective potential Ui{R) + 
U2{R) within this range of variation. Figure Ofb) shows this for the specific case of the 
Volkov-1 interaction. 

In summary, we have explored a model of the tetraneutron as a dineutron-dineutron 
molecule. Using a variational calculation we have found an effective Schrodinger equation 
for the relative motion of the dineutrons, after a proper account for the Pauli exclusion 
principle. An effective potential for the relative motion of the dineutron molecules was 
obtained. We showed that this potential does not have a pocket and thus the tetraneutron 
is very unlikely to be bound as a dineutron-dineutron molecule, although more complex 
variational approaches still can be explored. For example, one might consider a different 
spatial wavefunction for the dineutron system. We have used Gaussian for convenience. But 
the asymptotic form of the wavefunction might not be correct. An Yukawa form might be 
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more appropriate to achieve a lower energy for the system. Another improvement might 
be the use of more reahstic interactions, other than the Volkov interaction. Interactions 
including tensor parts might lead to some modifications of our results. Finally, although 
less probable as an improvement, one might relax the assumption of singlet states for the 
dineutron allowing for the triplet configurations in the calculation. 

Our study is complementary to other approaches and reinforces the commonly accepted 
idea that a tetraneutron is not a possible outcome of a theoretical calculation starting 
with underlying two-body nucleon-nucleon interactions. If the tetraneutron is bound, most 
probably it will be due to a special four-body attraction in T = 2 states or an exceptional fine 
tuning of the nucleon-nucleon interaction which does not seem to fit within present nuclear 
models. It might however be useful to recall an old phenomenological argument Q against 
the stability of the tetraneutron. Adding a pair of neutrons to a nucleus one usually increases 
the separation energy of the proton. If this rule holds, a simple comparison of the particle- 
stable tritium and unstable immediately leads to the conclusion that M^^n) > 4M„. 

After this paper was completed, a detailed variational Monte-Carlo calculation for the 
tetraneutron was done in Ref. It was shown that it does not seem possible to change 

modern nuclear Hamiltonians to bind a tetraneutron without destroying many other suc- 
cessful predictions of those Hamiltonians. Otherwise, our understanding of nuclear forces 
would have to be significantly changed. 

More elaborate QCD calculations based on lattice gauge theories cannot presently assess 
this problem. There are some lattice calculations for the H-dibarion state with not very 
reliable results (see Ref. 0, Q). However, lattice claculations still cannot determine if, 
e.g., the deuteron is bound. State of the art lattice numerical calculations are aiming at 
(9(10) MeV accuracy in binding energy, not 0{1) MeV or less that is relevant for nuclear 
physics. Also, how the calculations depend on quark masses is an interesting but very 
difficult question. 
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